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SUMMARY 


The present paper derives Von Karman’s Fourier integral method in 
supersonic wing theory directly from the basic concepts of the harmonic 
source and doublet. The method is first applied to investigate the 
general solution of the wave drag of a tapered swept wing with a sym- 
metrical diamond airfoil profile. The general solution includes all 
kinds of wing plan forms which may be swept backward or forward, and 
tapered or reversely tapered to any ratio. A number of the limiting 
cases are also investigated. For practical aerodynamic design, two 
families of wing plan forms with the fixed taper ratios 0.2 and 0.5, any 
swept angle, aspect ratio, and Mach number are shown in graphs. Some 
particular applications are illustrated. 

The reversed-flow theorem on wave drag as shown by Von Karman and 
Hayes . checks well with the consequence of the general solution. This 
method shows a certain elegance as no conical-flow assumption is needed 
and the mathematics is powerful enough to obtain a general solution 
covering all possible geometrical arrangements without detailed 
considerations. 

While in recent years the direct problem of finding the lift dis- 
tribution with a given angle of attack of the wing has been well solved 
by the method of conical flow and by other methods, the present treat- 
ment, on the other hand, investigates the inverse problem, that is, to 
find the downwash distribution in the plane of a wing with a preassigned 
pressure distribution. This is particularly favorable with the present 
method. • The general solution of the downwash of the tapered swept wings 
is derived for the case where a constant lift distribution on the -wing • 
is preassigned. Of course, the method may be applied to- any lift or 
pressure distribution along the wing chord and span. The corresponding 
angle of attack of the wing and the downwash can be determined everywhere 
in the plane of the wing. To demonstrate the downwash distribution as - 
given by the general solution, the downwash of a number of wings, 
including a sweptback tapered wing with supersonic trailing edge and a 
triangular wing, is shown in graphs. 
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The general solution of the direct problem in the present formula- 
tion is a singular integral equation of the Weiner-Hopf type and is 
difficult to obtain. 


INTRODUCTION 


In 1935 Von Karman first showed that the concept of the Fourier 
integral can be adopted to explain the similarity of Prandtl wing theory 
and the theory of planing surfaces (reference l) . In 1946, Von Karman 
introduced the Fourier integral method to the supersonic wing theory 
(an unpublished report of Northrop Aircraft, Inc.). The present author 
under Von Karman' s guidance investigated this problem to a certain extent; 
some of the results have been published (reference 2) but some, presented 
in 1948 at the Sixth International Congress for Applied Mechanics, are 
not available in published form. The present report may be considered 
as an extension of the earlier work. As some of the literature is 
unavailable at this moment, some of Von Karmen's work is repeated. 

As far as the linearized theory on supersonic wings is concerned, 
the conical-flow method was originated by Busemann (reference 3) and 
later reproduced and extended by Lagerstrom (reference 4), Hayes 
(reference 5), Stewart (reference 6), Laporte and Bartels (reference 7), 
and Snow (reference 8). Each of the later investigators approached the 
same method with slightly different techniques but confirmed the essential 
results of Busemann. Before the conical-flow method was known in this 
country, Jones visualized the advantages of the conical flow and showed 
some of the basic physical concepts (reference 9) in 1945. Later, in 
1946, with the concept of conical flow, Jones also showed the invariance 
of the Lorentz transformation and introduced the oblique coordinate in 
the swept -wing problem (reference 10). 

The challenge of supersonic flight has aroused the interest of many 
other investigators. Puckett showed that the source integration method 
could be applied to study the wave drag problem (reference ll) . Later 
Eward extended this method to solve the lift problem of a finite wing 
of any plan form (references 12 and 13) . Heaslet, Lomax, and Jones 
extended Volterra's and Hadamard's method to the supersonic wing problem 
(references 14 and 15). Gunn applied the operational calculus to the 
same problem (reference l6) . There are many other works, all of which 
cannot be mentioned in this paper. 

On the other hand, the basic concept of the Fourier integral method 
is quite different from the above methods. Instead of the concept of 
conical flow or simple sources and doublets, the present method considers 
along the direction of flight infinitely long harmonic source (or doublet) 
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lines, the behavior of which is quite equivalent to the harmonic acoustic 
source (or doublet) in the sequence of time. Each of such source lines 
will send out a divergent cylindrical wave in the radial direction. The 
potential of such a wave can be expressed by the product of the source- 
strength function and the Hankel function of the second kind. For each 
eigenvalue or frequency of the source oscillation, there is a corre- 
sponding eigenvalue in the argument of the Hankel function. By means 
of the principle of superposition, an arbitrary distribution can be 
synthesized with such simple harmonic sources of different frequencies, 
if the frequency spectrum is a continuous one. The most powerful tech- 
nique to serve such a purpose is the Fourier integral method. 

Now in the finite-wing problem, the boundary is considered as a 
source sheet rather than as source lines. Such a sheet can be built up 
by integrating elements of a source strip which is equivalent to a 
source line in behavior. 

The above simple physical interpretation may help to give some 
insight into the mathematical theory to be discussed later. 

The present paper first tries to give an introduction to Von Karman's 
method in terms of the physical rather than mathematical approach. Then, 
the wave drag of a tapered sweptback wing is treated. This work is not 
a duplication of that of some other concurrent investigators, but rather 
a supplementary contribution. The next step is the lift problem. For a 
given angle of attack, on an arbitrary wing plan form, it is rather 
difficult to solve the integral equation analytically. But for a given 
lift distribution the downwash angle can be evaluated anywhere in the 
plane of the wing.l 

The method has a certain elegance. The conqplete physical effect 
of the arbitrary wing can be expressed in terms of one integral. The 
finite number of discontinuous points of the Fourier-Bessel integral 
gives exactly the right picture. 

This investigation was conducted at the Department of Aeronautics 
of the Johns Hopkins University under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
Due appreciation should be given to Misses V. O'Brien and P. darken for 
their assistance in carrying out the research. 


-*-In some recent unpublished Navy reports, J. D. Miles has obtained 
a solution of the direct problem, particularly for the rectangular wing 
tip with constant angle of attack. 
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GENERAL THEORY OF THE FOURIER INTEGRAL METHOD 
Elementary Solutions 

In the three-dimensional, steady supersonic flow of a compressible 
nonviscous fluid, the differential equation of motion in the linearized 
sense is 

0yy + 0zz = (« 2 - l)0xx ^ 

where the velocity potential _0 is defined by 

grad 0 = q - Ui = ui + vj + wk, the disturbance velocity vector. It 
may be considered as a two-dimensional wave equation if x is conceived 
as time in the sequence of which the future can contribute nothing to 

x 

the present. This can be shown clearly by introducing t = —7—; 

(m 2 - l) 1 / 2 

that is, 

0yy + 0zz = 0tt (2) 

with the dimension in the flow direction being equivalent to time t. 

After this transformation, the velocity of propagation in equation ( 2 ) is 
unity. In this sense, equation ( 2 ) becomes the potential problem of the 
acoustic source and doublet in the y,z-plane in the sequence of time as 
pointed out by Von Karman (reference 2 ). After this transformation, the 
wing plan form has to be readjusted so that the Mach waves are inclined 
backward at 45° from the flow direction. As an example of transforma- 
tion, figure 1 gives a sweptback wing in the physical and transformed 
planes . 

✓ 

For a simple time-dependent harmonic source with strength cos Xt 
(where x/2n is the frequency of oscillation) located at (tj,0) in the 
y,z-plane, the elementary solution of equation (l) as given by Lamb 
(reference 17 > p. 297) is 

0 O ( t >y> Z J Tb°) = Re jj | e iX,t H</ 2 )(\r)J (3) 

where r = ^(y - q) 2 + z 2 and H,~^ 2 )(Xr) is the Hankel function of the 
second kind of zero order. It is easy to show that the solution satisfies 
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both the differential equation (2) and the 
cos Xt = Re e'*'^. At very large values of 

0 o (t,r) £ 1 cos[~X(t 

y 8it\r L 

the amplitude of which is inversely proportional to the square root 
of r. Physically, this represents the potential at the point (y,z) at 
time t due to a harmonic source of strength cos \t at distance r 
from it. Such a disturbance potential is called a divergent wave. This 
disturbance propagates from the point (t),0) with a circular or cylin- 
drical symmetry. 

Now if the strength of the source is the real part of F(\)e ±} ^ 
instead of just e*^, the potential must be also written as 


0 o (t,y,z', * 1 , 0 ) = Re i F(X)e iXt H 0 ( 2 ) (Xr) ( 4 ) 


given boundary condition 
Xr equation (3) becomes 



where 


F(X) = F q (X) + iF x (X) 


(5) 


with Fq(-X) = Fq(X) being even and F;l(-X) = -F^X) being odd. It 

is understood that F(X) is a function representing the strength, 
amplitude, and location of the source line. 


If two such harmonic line sources of equal but opposite strength - 
the negative one located at (q,-<;/ 2) and the positive one at (q,£/2) - 
approach each other, a harmonic doublet line can be obtained if 

lim F(X)£ = G(X) is considered as a finite quantity: 

- 

0l(t,y,zj q,0) = ^lM Q |0 o ^t,y,z; q,|) - 0 o (t,y,z; q, 


t- ^°r 

= iim — ^ 

i—>o 




= Re I - ^ GtXje^H^^Ur) 


5 ] 


( 6 ) 
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where the doublet is defined as positive, and 

G(X) = G 0 (X) + iG x (X) (7) 

with Gq(-X) = Gq(X) being even and G^(-X) = -G-^(X) being odd. It is 

very interesting to see that this potential of the doublet no longer has 
cylindrical symmetry, but is antisymmetrical with respect to the 
t,y-plane. The harmonic doublet line will be used for the wing lift 
problem. 


Boundary Conditions in the Supersonic Wing Problem 


Let a flat body or wing occupy a region in the x,y-plane and a + 
and a_ be the slopes of the upper and lower wing surfaces at (x,y), 
respectively. Both are small in comparison with unity and are zero 
outside the finite region occupied by the wing. Within the approxima- 
tion of the linearized theory, the wing may be considered equivalent to 
the superposition of a symmetrical body and a mean-cambered surface . 

The slope of the symmetrical body is 


a o = |( a + - a -) 

and the slope of the cambered surface is 


(8) 



(9) 


where a Q and a-^ are considered at the x,y-plane. In other words, 
the effect of exact location of a Q and a-^ in the z-direction is 

entirely neglected in the future treatment which has been shown by 
Von Karmen to be consistent with the linearized theory for the flat body. 
From the above, the required boundary condition of the potential is the 
normal velocity to the wing surface which is zero; that is, 

qn = 0 (lO) 


In this case of a flat body, where the vertical velocity component pre- 
dominates, equation (10) can be approximately expressed by 



= Ua + 


z=+0 


( 11 ) 
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w 



( 12 ) 


where w + and w. represent vertical velocity components on the top 

and oottom of wing surfaces, respectively. Besides, at a distance far . 
away from the wing in the upper stream, u=v=w =0 if there is no 
other disturbance generated in the upper stream. 


Supersonic Flow about a Flat Symmetrical Body 

If a source line at (tj, 0) is of any arbitrary strength with finite 
time t the potential can be built up with equation ( 4 ) by means of a 
Fourier integral as 



where U is the free-stream velocity. The insertion of U in the 
source -strength function is purely for future convenience. 

The above integral and its first derivatives exist if the Lebesgue 
integral 


Hoo 

|FU,Tj)Ho ( 2 )(\r)| dX 
UO 


< OO 


The above condition is automatically satisfied, if the source with finite 
strength is in action for a finite time interval. Actually, equation ( 13 ) 
can be considered as the potential of a body of revolution. Sometimes, 
in conforming with the convention of the complex Fourier integral the 
above equation may be conveniently written as 


0 o '(t,y,z) 



dXF(\,q)e m H 0 ( 2 )(Xr) 


(1M 
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with the understanding that the path of integration of X is slightly 
above or below the origin in the complex. X-plane so that the singularity 
at X = 0 can be detoured. It can be verified easily because 


F(-M) = F q (X,ti) - iFj(X, T]) = F(X,tj) 

Ho (2) (-Xr) = H 0 (2) (xre in ) = (xre' 1 *) = E^ 2) (\r) 


where the bar represents the complex conjugate of the original function. 
Now, the above source line may be considered as an elementary strip of 
a source sheet in the t,y-plane of width dT) with the above strength as 
the strength density per unit area. It is clear then that the. potential 
due to a source sheet built up by such strips is 


0 o U,y,z) 


^00 

0 O * d, l 




-00 



poo 


u 


— 00 


poo 


dt] 


dXF(X,T})e iXt H 0 ( 2 )(Xr) 


( 15 a) 


For the convenience of the later development, f(X,Tj) = - i F(X,t]) is 
introduced into equation ( 15 a), and then 


0o(t,y,z) 



dXf(X,n)e iXt H 0 ( 2 )(Xr) 


(15b) 


One thing should be pointed out. Owing to the structure of the integral, 
the source sheet may be composed of a number of discontinuous portions 
in the y-direction, and the source -strength function f(X,q) may be 
different in different discontinuous portions, but it must be the same 
within one continuous portion of the source sheet. Of course, in that 
case, the integration with respect to q has to be broken up accordingly. 
Physically, this integral can give the potential not only due to one wing 
but also to a number of wings or flat bodies in the t,y-plane. There- 
fore, it can be used to study the interaction between wings. 
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The preceding equation will he the disturbance potential of the wing, 
if the relation of the source-strength function Re f(X.,T|)e : '-^' can be 
identified with Oq or the thickness distribution of the symmetrical 

flat body. This can be obtained from the boundary condition. Differ- 
entiating equation ( 15 b) with respect to z, there results 


d0 

_o 

dz 


Ui 

2 



X dXf(\,q)e 1Xt H-.( 2 )(Xr) - 


( 16 ) 


Now as z — *0 the integrand goes to zero except in the neighborhood of 
q = y. Let q = y + z tan 0; dq = z sec ^6 d 6 , where 9 = tan - -*- ^ 

z 



y n 
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From equation (10), 


or 



From the theory of complex Fourier transform, 

J OO 

o-(t,y)e“ iX,t dt (19) 

00 

This means that, if the source distribution f(X,y) is chosen to be 
the Fourier transform of the distribution of the angle of attack, the 
potential 0 O will represent the flow about the wing automatically. 

Since the angle of attack is zero ahead of the wing, it can be shown 
easily that equation (15b) satisfies the second boundary condition of 
equation (12). 


Supersonic Flow about a Lifting Surface 

By means of the same argument as in the last section, the potential 

can be formulated for a doublet sheet of strength Re e^^G(X, tj) from 
equation (6) : 


0l(t,y,z) 


Re 


-iU 

T 



dX X G(X,Tj)e iXt H 1 ( 2 )(Xr)- 


( 20 ) 


where a constant, free-stream velocity U is introduced. In order to 
simplify the later development, further introduce 
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g(M) = = go(^il) - IgiUjT]) 


where, with equation (7), 


g 0 (M) = — ^ — (even) 
8m' 


XGr\ ( X , *n ) 

gl(^n) = (odd) 

8 m' 


Then equation (20) can he written 


'oo n oo 


0 x (t,y,z) = Re -2m'U / dq dXg(\ > q)e i ^ t H 1 ^ 2 ^ (X.r)^ 

U- 00 (Jo 


The above equation can be written entirely in the complex form 


>oo n oo 


0l(t,y,z) = -m'U I dT) dXg(X,Ti)e iX ' t H 1 ( 2 ^(X r )i 


' - oo U- oo 


if the path of X is chosen below X = 0 in the complex X-plane, 
because 


g(-^ii) = go(^ T i) + igi(^Ti) = g(^n) 

2 ) ( -Xr ) = H^^Xre -1 *) = ^(^(Xr) 


The path of X with Re X > 0 cannot he chosen in order to 
detour X = 0. In that case. 


( 2 )(-Xr) = Hi( 2 )(xre in ) = -Hj_( 2 )(Xr) - 2J 1 (Xr) ^ H-^^Xr) 


which is impossible to reduce back to the original equation (23). 
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In order to find out whether equation (24) for the doublet sheet 
satisfies all the boundary conditions of the wing or not, first differ- 
entiate it with respect to t. 


a0i 

St 


m 1 


dx 


= m'u = -m'Ui 




dXg(X,,Tj)e^^Hj^^ (\r)^ 


(26) 


The above equation is identical in form with equation (l6) if 
f(X,il) is replaced by ?L-g(X,Ti). Therefore, with the same technique, 
it can be shown that, as z — >+0, t] — > y. 



-m'U 

2 


dXe iXt g(\,y) 


(27) 


Now in the linearized theory the pressure coefficient on the wing 
surface is 


r (+ .r\\ -2u -2 

C p(t,y,+0) = — = — lr— I 

• U Urn' \dt / z= +o 


(28) 


Thus, with equation (27), 


Cp(t,y,+0) 


d\e iX ' t g(X,y) 


(29) 


With the Fourier transform. 


gU,y) - ± 


dte 1Xt C_(t,y,+0) 


(30) 


Therefore, if the pressure distribution along the wing chord is 
given, the doublet distribution function g(X,y) of frequency \/2n 
can be obtained with equation (30). The function g(X,y) is inde- 
pendent of the pressure distribution of the neighboring chord. Further- 
more, if A,— >0 in equation ( 30) > 
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g (°' y ) = / <itCp(t,y,+0) 


( 31 ) 


On the right side, dtC p (t,y,+0) is equal to one-half of the lift- 


-l/^ C L) 


coefficient distribution, — U — •), at y. Thus, g(0,y) is equivalent 
* v ' y ' 

at y which has a physical meaning. 


to 


Unfortunately, equation (23) of the doublet-sheet potential satisfies 
one of the boundary conditions on the wing, but does not satisfy the other 
required boundary condition that the potential and its derivatives must 
be zero far ahead of the body. In other words, the correct potential 
must be zero as t — »-». 


The potential of the doublet sheet at large ±t is investigated as 
follows. Introduce v = Xt into equation (23): 


0(t,y,z) = Re 


»00 'OO 

-a»'u an f 


For very large values of t there can be written 


Hi 


( 2 ) 


(t) 


2 it 

rtyr 


and 


~g(o,Tj) 


lim 0q(t,y, z) = lim Re 

t — } zb co t ^+00 


00 iv 




dv 


= lim 

t — > ±00 


Wd 


Jt | T c 
-00 


dTlg(0,Tj) 


sin V 


dv 


(32) 


0l( ±oo >y> z ) = X* = ±2m'U I dT]g( 0 ,T]) 


(33) 
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where the negative sign corresponds to t =•-<» and the positive sign 
corresponds to t = 00. As x (y , z) is independent of t, it satisfies 

-s2 * ^2 y* 

the Laplace equation - x + d --- = 0 and can be added to j, the 

<3x 2 By 2 

doublet-sheet potential, without an effect on the pressure distribution. 
It should be noted that the value of the potential function as t — >00 
is 2X*. Hence X* may be interpreted physically as the potential 
function of the downwash in the Trefftz plane. 

Thus, from equations (24) and (33) the complete solution of the 
present problem on supersonic wing lift is 


0(t,y,z) = 0i + X* = m'U 


\ dr] 2g(0,q) - r dXg(X,r 1 )e iXt H 1 (2 )(Xr) 

J-» r L J-co 


In the case where Cp is given, the above equation can be written, 
with equations (30) and (3l)j 

0(t,y,z) = f -%dr) | dTC p (T,i|,+0) (2 - f r dXe iX ’( t ‘ T *H 1 ( 2 4xr)] 


2 " J-.r 2 


Let Xr = v; then 


0(t,y,z) =|^ hq <irCp(T,q,+0) 2 - I dve 

(J- 00 J-00 U-00 


{¥) 




With reference l8, page 405, the bracketed term in equation (36) can be 
evaluated as 


V 2 ;, (v) - 


It - TJ 


t zJL < 1 


t - T 
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Thus, if the vortex of the foremost Mach cone is taken as the origin, 
there can be written 


0(t,y,z) 



(38) 


which is the equation of the potential in terms of the pressure coeffi- 
cient on the wing. The value of C p is zero outside the wing. Since 

the radical must be positive, 0 is zero outside the foremost Mach cone, 
as predicated by other theories. 


AERODYNAMIC BEHAVIOR OF SYMMETRICAL FLAT BODIES 
IN SUPERSONIC FLOW 


As the disturbance potential for the symmetrical flat body in 
supersonic flow has been given in equation (15b), with the Fourier 
integral method, the present section will show the expression for the 
pressure coefficient and the wave drag of various types of wing plan 
forms and airfoil sections. 


Derivation of the Expression for the Wave Drag of a Wing 

To the order of approximation of the linearized theory, the incre- 
ment of pressure Ap anywhere within the disturbance area is 


Ap = p 


80 o 

- Po = -pUu = -p Q U = 


PqU 

m' St 


(39) 


where the relations x = ( M 2 - l)V2 t = m't and u = 80 o /8x are used, 
p and p are the local pressure and density, respectively, at the 
point (x,y,±z), and p Q and p Q are the free-stream pressure and 
density, respectively. 

To find S0 o /dt, equation (15b) can be differentiated as 


u d0o 
m' St 



X dXf(X,T))e iXt H 0 ^ 2 ^(>.r) 



is an even function with respect to z. 


( 40 ) 


where 
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Substituting equation (40) into equation (39), 

Poo 

dq | X d\f(X,q)e iX ' t H 0 ^ 2 )(X.r) (4l) 


Ap(x,y,z) = 


PnU £ 


2m' 


If the pressure increment on the wing surface which is located 
equivalently at z = +0 is desired, r— »|y - q| and 


£p(t,y,±0) = 


PpU g 

2m' 


U 


Poo 

dq 

■ 00 J 


X dXf(X,q)e iX ' t Ho^ 2 \x|y - q| ) (42) 


If it is the wave drag distribution along the span cSD/6y which is 
of interest, all the horizontal components of the pressure increment 
along the- x- or t-direction on both the top and bottom surfaces over a 
uni t span width must be summed up. Mathematically, 


= 2 J"* Apa 0 dx ( 43 ) 

where a 0 is the surface slope of the airfoil along the line of con- 
stant y as given in equation (l8). With the expression for ao(t,y,+0) 
and £p(t,y,+0) in equation (42), the drag distribution per unit span 
can be determined as 


'00 f ’00 


p 0 u 2 dq X d\H 0 ^ 2 ) (X|y-q|)f(X,q) 

L/-00 vj - 00 



dX'e i(X+V)t 


f(V,y) 


(44) 


where X' and y correspond to Oq. To evaluate 
integral, write X* = -X" and the double integral 


the last double 
can be written as 
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dt I" dX'e i(X+V)t f(V,y) = | dt I (-dX")e i(X " X ’" )t f(-V , ,y) 


' -00 -00 


U -oc U -00 


dt dX."e 1 ^“ >L "^ t f(X.",y) 


where 


f(-X",y) = f 0 (-X",y) + if!(-X",y) 


= f 0 (X”,y) - ifi(X",y) 


= f(X",y) 


is introduced because fQ is even and fj is odd. From the complex 
Fourier integral theorem. 


2«f(\,y) 




dA." e i( X " X ")tf(x.",y) 


(46) 


if the Lebesgue integral 



dX." |f(X.",y) | 


< 00 


exists or if 


fU",y) 


is of bounded variation, in the neighborhood X = X" . (See theorem 23, 
reference 19, p. 42.) Substituting equation (46) into equation (44) 
gives 
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= 2np 0 U 2 J* dr\ J*° X " T| | ) f ( X, ti ) f( X, y ) 

poo poo 

= 2np 0 u 2 I dq I X dx|~Ho( 2 )(X|y - q| )f(X, q)f(X,y) + 

U-co U 


Ho( 2 )(x|y - tj| )f(X,q)f(X,y)] 


= 2«p, 


poo poo p 

Q U 2 dT) x dX J J 0 (X|y - r\\ ) [fo(X,q)f 0 (X,y) + 

J-00 Jo ^ 

f i( X, Tl ) f x ( )] + Y 0 (X|y - q|) [fi(X,q)f 0 (X,y) - fo(X,q)fi(X,y)[J 


(47) 


where the relations 


Hq( 2 ) ( -X|y - T| | ) = H 0 ( 2 )(x|y - nl) = Jo(^ly - nl) - iY 0 (X|y - T) |) 


f(-X,q)f(-X,y) = f(X,q)f(X,y) 

are used so that the integration limit changes from 0 to 00 instead 
of -00 to 00 . The imaginary part automatically cancels out. 


Now, if D, the wave drag of the whole wing system, is desired, 
equation (47) can be integrated with respect to y: 

^ dJ] I- X dXH o (2)(x,y " ^DfU^KU^y) 

= 4np 0 U 2 I dy J dq I X dX Jjo(X|y - ti| ) j^fo(X, T])fo(\,y) + 

J— 00 J-00 (JO ^ 


f 1 (X,n)f 1 (X,y)] +Y 0 (4|y - ill)[fi(X,Ti)f 0 (X,y) - f 0 (X,q)fi(X,y)]j. 


( 48 ) 
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Now, as the wave drag D is independent of y and Tj, it should be 
independent of the order of integration or of the interchange of y 
with T). In carrying out such an interchange, 


noo 


D' = 4jtp 0 U 2 


<*n 



- y|) [f 0 U,y)fo(M) 


fl( 4,y)fi(X, 



+ 


YoUh 


- y|) [fi(4,y)f 0 (X.,r|) 



This shows that the necessary condition D = D 1 is 



100 

ay 

-00 u 


poo 


HU 


Y 0 U|y - r il) |£i(^^y) f o( x ni) 


fl(^Tj)f 0 (X 


,y)]| 


= 0 


(49) 


This is one of the important relations that has been developed by 
Von Karman. 

There is another way of evaluating the wave drag as pointed out by 
Von K^rm^n. The total momentum transfer through a cylindrical surface 
of very large radius R which contains the wing in the flight direction 
must be equal to the wave drag, by means of Newton's laws. Thus, 


D = 




d 0 MM 

dr dt 


(50) 


which is the same as equation (48). Applying equation (4-9), equation (48) 
can be written as 


D = 4jtp 0 U 2 I dy I dT) I X dXJ 0 (X|y - q | ) Jf Q ( X,q)f Q ( X,y) + 

U -00 U -00 

fl(X,Tj)fi(X,y)J 


(51) 
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In the case of a sweptback wing which is symmetrical with respect 
to the x,z-plane, fo(X,y) and ^i(X,y) are even functions with 
respect to y. To take care of this point, write 


^ 00 poo poo 

D = W 0 U 2 dy dq X dXJo(x|y - q |) [f 0 ( X, \ q | )fo( X, | y | ) + 

—00 U-00 U 0 


fiU>hl) f i(MylO 


n 00 Poo Poo P 

= knp Q U 2 J dy J dr, J X dX^J 0 (X|y - q|) + 
J 0 [x(y + -n)] j [f 0 (M)fo(x, |y|) + fi(M)fi(My|)] 

r oo poo n 00 

dy j dr t j X dX-jjo(X|y - r,|) + 


JO JO 


J 0 [x(y + q)]| [f 0 (X,r,)f 0 (X,y) + f x ( X,q)fi( X,y)] 


(52) 


where the limits of r, and y change as shown. Equation (52) will be 
used in the next section very often. If the span is finite ( say that 
b is the semi span) and the airfoil sections are similar everywhere, 
equation (52) can be written 
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Wave Drag of a Tapered Sweptback Wing 

If the double-wedge profile is taken as the airfoil section of the 
wing as shown in appendix A , and equation (a 8) is introduced into 
equation (53 ) , there results 


Pb Pb 

D = 8jtp 0 U 2 dy dT] 

Jo c 0 



{ 


X dX •< J 0 ( X| y - T| | ) + 



cos 


Xay')(! 


cos Xa.-q * ) cos Xp(|y| - | T) | ) 


8 PqU 2 (%') 2 

n 





Jo [>( y + Tj )|J ( 1 - cos Xa y , )(l - cos Xa^cos X0(|y| - |q|) (5^) 

By definition. 


V .-V 

y' =■ 

v =4 

v* =^r 



a Q 

a o 

m ? 



ay f = a Q 1 ( 1 

- v'y') 

1 

> (55) 


a T1 t = a 0 f (1 

- vV) 






For details of the notation refer to figure 1 and appendix B. 
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The wave drag equation can be written in another form 


nD 


Pb' 


8p 0 U 2 (a 0 ') 2 (a 0 ’) 2 Jo Jo Jo 


Ob' 


dy' 


Hoo 


dtp 



Y ^Jo[(n’ - y’ )^o'J + 



Jo [(V + y* )Xa 0 £J |(- r ] 1 - cos ^o'jf + PU' - y') 


v 1 ( T) 1 + y' )J i jl - cos Xa 0 ’[2 - P( Ti' - y' ) - 
v’(V + y')J ^ - ^Jl - cos Xa 0 ’(p + v’)(q' - y* )J 


I jl - cos [> 


Xa'(p - V * ) ( T 1 * - y’fQ + 


i|l - cos Xao'jl + Pill' - y' ) - v'y'J^j- + 

y ») _ v'y’J^j- + 

y’ ) - v-n']j + 

3 }‘ 


1 jl - cos Xa 0 ’ [l - J3( T> ' 
||l - cos *v[l + P(ti* 
i jl - cos Xa 0 r |l. - P(ti' 

[l - cos p\a 0 '(V - y' Fj^ 


- y r ) - vV 


(56) 
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With the known infinite integrals of Bessel functions as shown in 
appendix C, the integration in equation ( 56 ) can he carried out. 


” D dy. dti'T- 1 cosh -1 2 ± Pin’ - y') - v'(n' + y') . 1 cosh -i a 4 p(y - y) - v -(ti' + y) 

8p 0 u 2 M 2 (ao’) 2 Jo Jo L 1* - y ' 4 I' ty ' 

1 cosh -i l - p(v-y) - v'(v*y) . 1 cosh -i g - my - y) z_hlLnI i y,) . 
k v - y' ^ V + y* 


4 


— cosh" 


cosh" 


cosh- 


V + y* 

4 

1 + 3(ti' - y’) 

- v*y’ 

V - y’ 


1 - PC 11 * - y’) 

- v’y’ 

n* - y* 


1 - PU* - y’) 

- v’V 


I (e -y')(V - y') 

n r + y r 


n' + y’ 


cosh" 


n* + y r 


n r - y’ 


n r - y* 


V + y’ 


V + y* 


n* + y' 

cosh-i p I (57) 


h- 1 pj 


In all cosir 1 (A) terms , A must be positive. In order to carry out 
the integration, a transformation of the coordinate system may he intro- 
duced as follows: Let u r = q 1 + y T , v ! = rj f - y T . Then, 


dy r dr] f 


S(y* ,T)' ) 

d(u f ) 


du f 


dv 1 



dv f 


= 0, 

V 

= o — 

u' 

= 0, 

v' = 0 

0 

II 

V 

= b — » u' 

= b. 

<; 

11 

= b. 

V 

= 0 - 

U* 

= b. 

v* = -b 

= b, 

ii’ 

= b — 

-> u’ 

= 2b' 

, v = 0 


With this transformation, the square domain in the y 1 ,^ 1 -plane will 
transform to a diamond- shaped domain (of twice the area) in the 
u 1 ,v* -plane. 
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_(b') 2 |jL cosh" 1 (0 ♦ v') + J cosh" 1 (0 - v') + cosh" 1 0 J.+ ^ du' ^ dv '^ 5 cosh" 1 ( 2 J u * 0) 
I cosh" 1 (i^l - v) - 1 cosh- 1 " 0) " i cosh" 1 (2^1 - »•) - t c^- 1 - 

1 cosh -i [- 1 cosh- 1 ( *-*■;; ♦ >)-* cosh- 1 (|^ - ,•) - 

t — ( 2 - 2 v T - - ") - i “-- 1 - 0-5 ““"- 1 - 

i — i f^-] * I [- 1 + ^ t) + -*- 1 [ 2t(2 ^ v ' )Y ' - t] * 


■1 1 

2 - v'u' O V'\ , h -l f2 - 

(23 - v' )▼' v'™| 

+ cosh- 1 I 2 • V ‘ U ' + P - v ') + 

V 

2v- P 2/ L 

2u* 2j 

V 2v’ 2 / 

- 1 r 

2 + ( 2 p - V* )v* v’l h -l / 

2 - V'u* v’\ 

. cosh- 1 P - (20 + W)V . v]l + 

L 

2u' 2j V 

2v' P 2/ 

L 2u* 2J r 


1 j— 1 ♦ .♦*■)* — 1 [ a ^:::r - t] + — 4 ( : 


2 - 2v , b l + y'u 1 
2v' 


[2 - (23 - v' )▼’ v'l 

_^_1 (2 -.2v'b’ + v’ 

'u* 

+ j 3 


. [2 + (23’ 

- v' )v' v’H 

L 2(2b* “■ u* ) 2J 

" \ 2v' 


2 J 

L 2(2b’ 

-u') 2 J 

(2 - 2v’b' + v’u* 

f\ + cosh-1 P - (20 + 

V* )v* 

-4 

1 + ( -cosh-1 — - 

cosh-1 <*' 

l 2v' ' P 

77 c0 h L 2(2b- - 

U') 

2] 

1 1 

2b' - u' 


The details of the integration are too complicated to he given here. 
The result can he summarized in equation (59) which follows, where 

3 0 = (3 + v' , Pp = 3 - v’, a = 1 - v'b' , and sin -1 (A) = (sign A)| 

if | A I >1. The upper line with 1 - 3^ should he used, if 3 < 1; 
the lower line with 3 2 - 1 should he used, if 3 > 1; and the signs 
of all cosh - -*- terms should he reversed if 3 < -1- The above rules 
should be also applied for 3 0 and 3-^ Therefore equation (59) covers 

all cases except a few limiting cases: 3 = 3 0 = Pi = 

v' = 0, and 3 = 00 (i.e., M — y 1) . These limiting cases will he shown 

in the next section. 
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* b’(i ♦ o) p- - (b 1 ) 2 

2 Cfc 


{^T7 I75 ““' 1 I^H] 'brrp^' 1 I^Q Ml 


a -»') z r . . k - t ' i 
|»(., 2 - i)‘*t I 1 - '»■ I 



. (» *«•..•) 2 L., 

«0 O » 2b’ 

, [..(■• 0 - 1 ] 

[“” * !ll • Jj 

2(fc •!)(>- «. 2 ) 1/2 l 

a ♦ 2B 0 b* 

L° • * 2 Bob’ Jj 

.L.fr-O‘1 

il. ("•»..■) 2 r 

th-l | aB o* ZV | 

+ I + (1 - R 

p »..»■] 

U 2( fc .,)( fc 2 -.)^l 

|o ♦ 20^* I 


r (. . j 

. L. 

o9 x - 2b 1 

.1 [ . ( l - pl2 M 

[ 2(1 - pv) J 

2 (». ■ >)fr • ». 2 ) l/2 1 

a - 2&ib' * " 

t 1 a - 29b 1 J 


2(9o 


1 1, (1 ’ ’ l2) ‘' Il 


[oPj - 2b 1 | 

L (‘:«i 2 )‘-|] 

| - co.h | 9l 2(1 . pb ,,|j 

2(Pi * 9)(&! a - l) V2 l 

|o - ZPib’J 

P 0 - 2Bb’ |J 

- 2v .1 r. (1 - p 2 )b'i 

| (# . 2Sb') 2 

Is Id ' 1 - slo’ 1 fp . 

J 0 * 29b’ [ 

{1 - 92) V l] - 

"ipV * l ° [ 9 0 - 2 P b' J 

J a(f>i ♦»)(! * p 2 ) 1 ' 2 

a * 20b’ j I 


O0 - 2b 1 I . c(jBh .i 

„ (l ' S 2 )^ 1 ’ 

t (0 ♦ 20b 1 ) 2 

rnflh -l \<# + 2V ! ^. h .! Ip . (l - P 2 )_b' 

[0 - 20b' 1 | 

| a - 20b’ 

2(0! > 0)(02 - l) 1 / 2 

|o ♦ 20b' j | 0 ♦ 20b 1 


_S_^ jv .... [„ - 4 [* - [«••(.- >•%] ■ fr- • fr - - if£b? ■•] 

^ 75 (f k • (1 ' y )b ' I - «* “•“• 1 k - I * T K * f 1 - ^7 1 - irh ““- 1 I* • fr - *?)'• | * ifoViy ““■ 1 w] 
[, ... [ fc . li^k] - , [, • . * - [, - c - ..1*] - ^ - [, - c - ,>] • ^ - *] 

2t , ( 6i a 1 _ i)i^ [-° g “"‘• 1 k • ^ ' S }t 1 * ^ c ° ,t ~ 1 k- (l 'V )l I ■ 7 1 ““‘ 1 h -fr 1** - - »!»' * ’)’ ”“‘ 1 n] 


{, ... [, . <^] . , .... [. - - fr • fr - A-] - a - fr.fr- -*•] • ’} 

i^rr^k 2 -* I* - • * — 1 1- • ^1 - irh — !• • fr - rtb 'l • irb -* I' • fr • ' 2) " 1 - (.“Vk) "“- 1 l# j 


where a, b, c, A, B, C, D, E, and F designate the terms to be 
used in the next sectipn. 

For a given sweptback wing, P 0 > 1, = 1, or <1 means the Mach 

line ahead, on, or behind the leading edge, respectively; 0j_ > 1, = 1, 

or < 1 means the Mach line ahead, on, or behind the trailing edge. With 
a fixed taper ratio, the trailing edge may be swept forward; P]_ = -1 
means the sweptforward trailing edge coincides with the Mach line. 
Similarly, p refers to the midchord or the line of the maximum thick- 
ness. At all p's equal to |l| there occurs a discontinuity of the 
slope of the curve of the wave drag. In other words, curves for P = 1 1 1 
are the envelopes of cusps of the C D /C Do curves at constant p. 
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To demonstrate the use of these results,, the wave drag coefficients 
of two families of sweptback wings have been calculated as shown in 
figures 2 and 3 with taper ratios of 0.2 and 0.5, respectively. The 
essential parameters are C D /C Do , A tan A, and 3 . The effect of the 
Mach number is contained in 3 = tan A/m’ and C Eo = ka^j m' , the two- 

dimensional wave drag coefficient of Ackeret. To use these curves , first 
fix the sweptback angle A and aspect ratio. Of course, a will have 
to be chosen as 0.2 or 0.5 in order to use these graphs. Then, with a 
fixed A tan A, read off C E /c Eo at various values of 3 along the fixed 

abscissa. Replotting C E /c Eo against 3, the example shown in figure k, 
actually gives a family of swept wings of fixed A tan A. The curve of C D 


against Mach number can be plotted from the relation 
This is not given in this report. 



It is very interesting to note that ( c D/ c Do)max occurs vhen the 
Mach line coincides with the leading edge, 3o = 1, or the line of maxi- 
mum thickness, 3=1. In figure 2 (a = 0.2), ( c D/ C Do) max occurs at 

3=1 for A tan A > 7-5 approximately, and at 3 0 = 1 for 
A tan A < 7 . 5 . In figure 3 ( cr = 0.5), ( c p/ c Do) may occurs always 

at 3o = 1 except for A tan A < 2 / 3 . Also, C E /C Do increases mono- 
tonically with A tan A at 3 0 = 1 , 3 = 1, or 3 -l = 1 , but decreases 
monotonically with A tan A for 3i = -1. 


Figure 5 shows 3 corresponding to 3o = 1 an( l Pi = 1 plotted 
against A tan A. This is needed for locating the cusp of the curve 
of C D /C Do against 3, as shown in figure 4. These curves are obtained 
from the relations: 


Pi = 3 


1 - 


2(1 - g) 

(l + a)A tan A 


1 + 


2(1 - a) 


(1 + a) (A tan a) 


l 


Po = 3 
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A few interesting facts can be drawn from figures 2 and although 
the taper ratios are fixed at 0.2 and 0.5j respectively. First, if 
p = tan A/m 1 = 0, M — ^00 for A / 0, or A = 0 for M > 1. For A / 0, 

Cp/Cpo = 1 or C D = c Do as indicated. For A = 0, 0 < Cp^Cp 0 ^ 1 

and the exact value of Cp/ Cj) Q is indeterminate. As p increases 
from 0 to 1, Cp/Cpo r i ses quickly with A tan A until it meets the 

curves p Q = 1 and then varies with A tan A very slowly. At p = 1, 
CpjCpo i s extremely high. This means that the closer the Mach line to 

the maximum-thickness line, the higher the drag. As p increases 
further, from 1 — > 00 , Cp/Cp 0 decreases with increasing p. This means 

that the closer the Mach line to the maximum-thickness line, the higher 
the drag also. 

If the sweptback angle is negative, or if the flow direction is 

reversed, the curve of Cp^Cp Q against -A tan | A | is actually the 

reflection of the present curves with the following changes: p — > -p, 

p 0 = 1 — } Pi = -I* Pi = 1 — ) 3 0 = “1; and - 3l = “1 — > 3o = 1 Po 

is defined as corresponding to the leading edge, and p^ as corre- 
sponding to the trailing edge. It will be found that Cp|Cp Q depends 

only on A tan | A | and |p| for the same wing whether it flies forward 
or backward. This confirms the reversed-flow- theorem in references 2 
and 5. 


Special Cases of the Wave Drag Equation 

Most of the special cases of the wave drag equation are required in 
calculating the curves shown in figures 2 and 3 and therefore must be 
written down one by one. These special cases are based upon the taper 
ratio that is between 0 and 1. 

(a) As Po — *1, expressions A and D in equation (58) take on the 
following limiting values while the rest of the equation stays the same 
except a = b = c = 0. 
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A = 


(1 - v'b’ + b' ) 2 

2(v' - 1) 

(1 - v'b 1 + 2b' ) 2 
2(v’ - 2) 


r 1 _ v'b' 

a - v'b 1 + b' 


/I - v'b’ - b'X 1 / 2 

\1 - v'b’ + b ’/ 

/l - v'b' - 2 b'\l /2 / 1 _ v'b' ' 

y. - v’b' + 2b ' j \1 - v’b’ + 2b'y 


D = 


2 v' 


a 3/ ^ 2 [(l - v’b' + b’ i 1 ! 2 - (1 - v'b' + 2b' ) ly/ ^j 


,( 60 ) 


2v' - 1 
v' - 1 


(1 - b' ) 1 / 2 - — — — — (1 - 2b') 1 / 2 + 


v' - 2 


(v’ - l)(v' - 2) 


This case occurs when the right characteristic coincides with the 
leading edge on the right side. 

(b) Similarly, for 3^ = 1, expressions B and E become 


B = 


(1 - v'b' - b’ ) 2 

/l - v’b' + b’ N 1 / 2 / 1 - v'b' N 1 / 2 

2(v' + 1) 

\l _ v'b' - b' / \1 - v'b' - b’ / 

(1 - v'b’ - 2b' ) 2 

(l - v'b' + 2b'\ 1 / 2 ( 1 - v'b' ^ 

2(v' + 2) 

\^1 - v'b' - 2b' y \l - v'b' - 2b ' j 


1/2 


-(61) 


E = 


2 v' 




1 - v'b’ - b 


)V 2 + (1 - 


V'b' - 2b 


) 1 / 2 ] . 


— + 1 (1 + b' ) 1//2 + 2 ^ V> + ^ (1 + 2b’ ) 1 / 2 + — 1 

v' + 1 v' + 2 (v' + l)(v' + 2) 


Here b = 0. All the other terms remain the same as shown in equa- 
tion ( 59 ). This case occurs when the right characteristic coincides 
with the trailing edge on the right side. 
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(c) When 3 ■= 1, the expressions C and F 


are 


C = 


(g + 2b' ) 2 

2(2 - v' ) 
(a - 2b' ) 2 


2(2 + v' ) 


~_ fo - 2b^W 2 + / a U /g] 

~\g + 2b'/ + \a + 2b'/ J 

/o_J_2V_\ 1 / 2 + / g X 1 / 2 ] 

\g - 2b'/ \g - 2b'/ J 


F = -ij- < g 3 / 2 [^-(g + 2b’ J 1 / 2 + (g - 2b' J 1 /^] 


r (62) 


2(1 - V' ) 

2 - v* 


(1 - 2b' )^ 2 + 2 ^ 1 + v ’ l (1 + 2b' J 1 / 2 - 

2 + v ' 


4v 1 


(J. - V 2 ) 


Here b = c = 0, All the rest remain the same as shown in equation (59), 
(d) When P Q = - 1* A and D become 


A = B in equation (6l), where 3^ = 1 


D = E in equation (6l) 


(63) 


Here a = 0. The rest of the terms in equation (59) remain the same. 
This case occurs when the left characteristic coincides with the leading 
edge on the right side. This conforms with the reversed-flow theorem, 

3 0 = -1 — ► 3p = 1, if the flow direction is reversed. 

(e) When 3p = -1, B and E become 


B = A in equation (60), where 
E = D in equation (60) 



(64) 


Here a = b = c = 0. The rest of the terms remain the same as given in 
equation (59)- 
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(f) When 3 = -1, C and F are the same as given above for 
3=1 with some rearrangement of terms and a = c = 0. 

(g) When v' = 0, or there is no taper, 3 0 = 3p = 3 and equa- 
tion (59) becomes 



The above equation is in accord with the results given in reference 2. 
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(h) If M — * 1, p — > ±<» as A is positive or negative. By 
definition, p = b/a 0 , v’b' = vp, = p tan A, <jd 0 = p tan A 0 , and 
= p tan Ap. 


- log log 2 + 

2 


56 • 0*b « (Hr) - - (H^)] * Hhh * - 

h • I - fl b Hf - « l‘ - HHrl) • 



There are a number of other limiting cases such as p o + p = 0 
(Aq = -A) and Pp + P = 0 (Ap = -A) which are also required in cal- 
culation. If a = 0 and Pp = 0, equation (59) will give the wave 
drag of the triangular wing which checks with Puckett's results. Owing 
to limited space, they are all omitted. 
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AERODYNAMIC BEHAVIOR OF A WING WITH A 
GIVEN LIFT DISTRIBUTION 


From the section Supersonic Flow about a Lifting Surface, 


gU>y) 



,-iXt 


Cp(t,y,+0) 


(30) 


Thus, when the pressure-coefficient distribution C p (t,y,+0) is given 
on the wing plan form, g(\,y) can be obtained with this equation. 

Consider a tapered sweptback wing with a constant lift distribution 


Cp(t,y) - Cp 0 
-b ^ y = b 


p 0 M - a o' = t = Pilyl a o’ 

where |3 0 and (Bp are defined in figure 1. . Thus, at any y within 
the span , 


g( X,y) 


^Pllylao' 

_1_ 

2 * 

tlpol y|-a 0 ’ 


-iXt ■ 
C po e dt 


■'■^po ~i^(Pi| y | a o' ) 

2*X L 6 


-ix(p 0 |y|- a o 

e 


i 


(67) 
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As X. — > 0, g(0,y) can be evaluated as 


c PO I , p o " P 1 , 


g(0,y) = — lao' -|y I 

n \ 2 


Substituting the above equation in equation (2h) , there results 


“ (t ' y ’ + 0) ■ T 2 J 


¥(v - - 


jy - nl f 00 dX f e - iX (Pi|y| +a o’) 
2 I 2«Xi L 

U -00 


s -ix(3„|y|-a 0 ')] e ixt Ei ( 2 ) (My _ M) 


-m* UC, 


(y - n)‘ 


Po - Pi, 


j ly^- nl ^ ^^“^^(Pi| y| +a 0 ') _ 

* 

e -lx(So|yl-ao')] e ^ Hl (2)( My . n ,) 
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and 


With the aid of appendix C, the infinite integrals can be evaluated 


v(t,y,+0) = 


- c po m ’ u an 


2 it 


(y - 


/ Po - Pi , \ 

\°' ' — 7 — hl j - 


|y-nl 


t i - Pi I’ll 

|y - nl 


+ 0, t r - Pilnl 5 |y - r|| 


*i - Pihl . 2 \( t i - ^l) 2 - (y ‘ , 

F^i ’ 41 • Bl '” 1 8 |y 


|y - q| 


to - Pohl 


|y - nl 


•> t 0 - pohl ^ |y - nl 


(70) 


tp - Pohl 2 \|(t 0 - Pohl) 2 - (y - 1) 2 


|y - n| 


|y - n| 


> t 0 - Poh I ^ |y - 'nl 


where t Q = t + a 0 ', t]_ = t - a Q ' . Only one of the two expressions 

holds for the range specified in the bracket. 
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If the integration of T| is written from 0 to b only, 


w (t,y,+0) = 


1 

°B 

U 

b 

d.T) 4 

2a o’ 

" (P< 

- Pi) 1 ! 

2a 0 ’ " (P 0 

- Pi) 7 ! 

4jt 

"J 

0 


(y - 

n) 2 

(y + t]) 2 

2a 0 ' 

_ (Po 

~ P l)^ , 

2a 0 ' 

- (Po - 

Pi)n 


\ 

1 

(y 

- n) 2 



(y + ti) 

2 


2 V(‘c 

) “ 

M) 2 ■ 

(y 

-T)) 2 

2 \fFo 

- Po 1 !) 2 - (y 

+ n ) 2 


(y - 

^) 2 



(y + t)) 2 


2 \[Fi 

- 

M) 2 ' 

(y 

- h) 2 

+ 

- Pi 1 !) 2 - (y 

+ TJ ) 2 


(y - 

Tl) 2 




(y + n) 2 

- 



"»b 

r «- 






- c po m ' 

u 


\l(to - 

• Pon) £ 

■ - (y - ■ 

T? 


2« 





(y - 

i) 2 




c 

0 






\l( t 1 

- sm) 2 - 

(y - 

n ) 2 

+ 


Po 1 !) 2 - (y + 

n) 2 


(y - 

Tl) 2 




(y + n) 2 



J] 


\ (^i " ^i^) 2 ■ (y + ti) 2 
(y + t}) 2 


( 71 ) 


where t 0 - ti - 2a 0 


is used so that the first terms cancel out. 
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Now, carrying out the integration with respect to t}. 


w(t,y,+0) 


Cpom’U 

2it 


[-(A - B) + (C - 



where 


pb 


A = 


dTj 


WO 


\|(to - (3 qti) 2 - (y - ~tQ^ 
(y - n) 2 


(see table I) 


(72) 


(73) 


Pb 


B = 


d-Ti 


uo 


\|(t x - Pxn) 2 - (y - 

(y - *i) 2 


(see footnote, table I) 


(7*0 


Pb 



\|(t 0 - M ) 2 - (y + n ) 2 

dn J 

(y + n) 2 


( see table II ) 


(75) 



^(t]_ - p^) 2 - (y + q) 2 
(y + q) 2 


(see footnote, table II ) 


(76) 


The present treatment is equivalent to considering two hypothetical 
wings. The first one is of constant positive lift distribution (C po is 
an assigned constant) starting at the leading edge and extending down- 
stream to infinity with semi span b. The second wing is of constant 
negative lift distribution (C po is an assigned negative constant) 



37 


NACA TN 2317 


starting at the trailing edge and also extending downstream to infinity 
with semispan b. The superposition of both hypothetical wings gives 
the aerodynamical behavior of the actual wing. 

Integral A, equation (73)> gives the contribution of downwash due 
to the right-side area of the first hypothetical wing in the forward 
Mach cone of the point P(t,y). Similarly, integral C, equation (75), 
gives the contribution for the area on the left side. 

On the other hand, integral B, equation (7*0, gives the contri- 
bution of the right side of the second hypothetical wing. Similarly, 
integral D, equation ( 76 ), gives the left-side contribution. 

As far as the limits of the integrals are concerned, each holds 
for within the limit 0 to b, whenever the integrand is a positive 
real quantity. However, for many locations of P(t,y), the integrand 
is positive only in a much narrower range of t) than the interval (0,b). 
In order to determine the valid range of the integrals and their respec- 
tive values for all possible locations of P, tables I and II are given. 
The corresponding values of the integrals for the different cases are 
given in the right-hand columns of the tables. 


Take table I as an example. If the leading edge is subsonic 
( 3o > l ) , the lower limit of T| is always zero because the forward 

Mach cone at P always intersects the wing center line within the wing 
area. When > 1 the upper limit of T| is given in the upper half 
of table I. Condition 1 concerns the location of the point P, which may 
be ahead of, on, or behind the leading edge, while condition 2 concerns 
the right intersection of the forward Mach cone of P with the leading 
edge, or side edge. Under both conditions 1 and 2, there are six 
possible cases. The upper limit of q is given in each case, and the 
value of the integral A is also given accordingly. 


If the leading edge is supersonic (3 0 < l), the integral A exists 
only when the point P is behind the leading edge. In particular, 
both the upper and the lower limit may vary within the interval (0,b). 
The condition for the lower limit depends on the intersection of the 
left side of the forward Mach cone of P. It is zero when and only 
when the intersection is on the center line of the wing in the wing 
area. If the intersection is on the right-side leading edge, the lower 
y - t Q 

limit is — as indicated. The condition for the upper limit 


1 - 3 0 

depends on the right-hand intersection of the forward Mach cone of P. 
If the right side of the forward Mach cone cuts the leading edge within 


the span, the upper limit is 
is b. 


tp + y 

1 + 3 o ; 


otherwise, the upper limit of T| 
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The tables for the other three integrals can be explained similarly. 
Figure 6 shows four typical cases concerning the integral limits. Fig- 
ure 6(a) shows the forward Mach cone of a point P(t,y) that lies only 
in the right-side area of the wing. Therefore only integral A exists 
and all the other three integrals are zero. Since the leading edge is 

/ \ t G + y 

supersonic (p n < 1), the upper integral limit is . It belongs 

V 1 + 3o 

to integral A-9 as indicated, according to table I. Figure 6(b) shows 
the forward Mach cone covering both the left- and right-side areas of 
the wing. Thus, only integrals A and C exist. The integrals are indi- 
cated. Figures 6(c) and 6(d) can be similarly explained. 


Although the lift distribution of the wing with supersonic trailing 
edge is well-known, the downwash may be interesting to explore. Fig- 
ure 7 shows three infinite half -wings with trailing edges at 15°, 30°, 
and 1+5° from the leading edge which is normal to the direction of flight. 
This is shown in the t,y-plane, or the x,y-plane with a Mach number equal 
to f2. A negative infinite downwash always occurs at the tip. The 
curves are plotted as 2rtw(y/t)/m' C p0 U against y/t. As is known, 

Cp 0 = 2a]_/m' , where is the angle of attack of the wing; therefore, 

nw/a^U can be plotted against y/t. Between the tip cone and the wing 

the downwash is constant but increases with increasing trailing-edge 
angle. Owing to the conical flow, the downwash is identically the same 
along the radial lines expressed in the conical coordinate y/t. The 
above curves were calculated with equation (72) and the tables for 
integrals A and B by setting p 0 = 0; Pj = tan 15°, tan 30°, and tan 45°; 
and a Q ' = 0. (Here t = t 0 = t-j_. ) 

Figure 8 shows a wing tip of unit chord with a raked angle of 30° 
for three values of t. At t = 0.5> it behaves exactly the same as 
the 30° case of figure 7(b), but at t = 1.5> it is quite different. 

Two infinite downwashes occur, one at the center of the leading-edge 
tip cone, and the other at the center of the trailing-edge tip cone. 

The former remains negative, the latter positive. Outside the two tip 
Mach cones, the downwash is identically zero as predicted by the two- 
dimensional theory. The above calculations can be made from equation (72) 
and the tables for integrals A, B, and C by letting f3 0 = Pi = tan 30° 
and a Q ' = 0.5* (Here t = t 0 - 0.5 = + 0.5- ) 


As another interesting example, the downwash distribution near the 
nose of an infinite constant -chord sweptback wing is calculated (fig. 9). 
The leading edge is swept back at 30° but is still supersonic. The 
downwash at different locations on the span is plotted along t. In 
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this case, the curves are plotted with 


at different values of y where Cp 0 = 


At y = 0, the downwash angle reaches the maximum value (0.8l65) 
at a single nose point and remains zero along the chord up to the 
trailing edge. There the downwash angle is negative and equal to the 
maximum value in magnitude and recovers to zero downwash downstream. 

For instance, at y = 0.4, the downwash angle is constant in the super- 
sonic region, and drops down under the influence of the nose Mach cone. 
At the trailing edge, the downwash angle drops to a negative value 
abruptly and continues to decrease until the trailing Mach cone is 
reached. The drop of downwash at the trailing edge is exactly equal to 
the downwash angle at the supersonic leading edge. Further along t, 
the downwash angle rises again and becomes asymptotic to 0 as t — >■ 00 . 
The same can be applied at any value of y until the nose Mach cone is 
off the trailing edge. 

Of course, it is very difficult to build a wing with the angle of 
attack as required by the constant pressure distribution. But by the 
present method it is so simple to calculate the downwash distribution 
anywhere in the plane of the wing that it will surely afford many 
applications if the trailing edge is supersonic. With a given pressure 
distribution, it is easy to calculate the sidewash; however, the details 
are omitted here. 

As another example, the downwash of a finite triangular wing with 
a 60° sweptback angle is calculated. Figure 10 shows the downwash dis- 
tributions at different values of y. Figure 11 shows the downwash 
distribution of a tapered sweptback wing. The taper ratio is l/6. 

Figure 12 shows the downwash of the same sweptback wing far downstream. 

It is very interesting to examine the effect of different Mach cones 
on the downwash. One important feature of the subsonic leading edge is 
the infinite downwash angle. As before, the above curves are calculated 
with equation (72) and the tables for integrals A, B, C, and D. 

The present approach to calculating downwash of a wing with constant 
lift distribution cannot be applied to wings with subsonic trailing 
edges, because the Joukowski-Kutta condition must be satisfied at such 
trailing edges. In order to investigate the case with a subsonic 
trailing edge, it is necessary to assume the distribution of the pressure 
coefficient drops to zero at the trailing edge. 


V 1 ~ ftp >f(t,y) 


a. 


2cu 


U 


against t 


m’\|l - 3o 2 


is introduced. 
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The preceding discussion just gives some cases with simple pres- 
sure distribution along the chord to demonstrate the application of 
the method. Of course the method is not limited to these cases alone. 


The Johns Hopkins University 

Baltimore, Md. , November 4, 194-9 
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APPENDIX A 


FOURIER TRANSFORM OF THE SLOPE DISTRIBUTION ON AN AIRFOIL 


As an example of the slope distribution on an airfoil, the Fourier 
transform of the slope on a double-wedge airfoil is calculated. 


z 



The accompanying figure shows the shape and slope of a double-wedge 
airfoil on a swept back wing at a distance y from the center of the 
wing span. The slope distribution is as follows: 


a(t,y) = 0, t <t D (y) - a'(y) 

a(t,y) = a 0 '(y), t 0 (y) - a’(y) < t ^ t Q 

a(t,y) = -ao'(y), t Q < t <t Q (y) + a’(y) 


(Al) 


a.(t,y) = 0, t >t Q (y) + a'(y) 
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where a, a 0 ', t 0 and a' are functions of y in general. Sub- 
stituting equation (Al) into equation (19), it yields 


f(x,y) = f 0 (x,y) + ifp( x,y) 


■fcj : 


d(t,y)e" iXt dt 


V e' 1Xt dt - V 


> t 0 +a' 


2 « 


lt 0 -a' 


2n 


-iXt .. 
e dt 


a, 


no 


e -iXt o| 
2n 


e~ iXr dr 


-a' 


Ha' 


JO 


e" iXT d T 


(A2) 


where t = t Q + t is used in substitution. It is easy to show that 


f(X,y) = — — e i ^° (1 - cos Xa' ) 
inX 


(A3) 


fo(X,y) = — sin Xt 0 (l - cos Xa' ) (even with respect to X) (a4) 


jiX 




f-,( X,y) = -2— cos Xt Q (l - cos Xa' ) 
x nX 


(odd with respect to X^ (A5) 
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For a rectangular wing with constant airfoil section, where a 0 ' , 
t Q , and a* are independent of y in the span, equations (a 4) and (A5) 
reduce to 


= 0 


f x ( ^y) 


a* 

( 1 - cos Xa 1 ) 

jtX 


if t Q is chosen equal to zero. 

For the case of a straight tapered swepthack wing, as shown in 
figure 1, 


a' (y) 


a o - y|y| 

m’ 


t o(y) = P | yj 


a 0 '(y) = a 0 ' independent of y 


For the details of the definitions of the different parameters, refer 
to the figure. Then, 


fo^y) 

f^U^y) 


sin (Xp|y|)(l - cos Xa y ') 

a n ' . . 

cos (X0|y|)(l - cos lay') 


(A 6) 
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Now, 


f(X,y)f(M) = 


-a 0 ' -iA.pl y | 

e 

ijtX 


(l - cos Xa y ' 




= e - 1XP(ly| - |, ' l) (l - cos Xa y -)(l - cos »«,') 

ji 


t a A 2 _ 

= — — 2 — ^1 _ cos Xa y ' H 1 - cos Xa^') Icos Xp(|y| - hi) - 

7i 


i sin 4 p(|y I - hi)] 


= [f 0 ( 4 ,y)fo( 4 ,n) + fi( 4 >y)fi(M)] + ijf 1 ( X,y)f 0 ( X,t) ) 


fo( x -^y) f i(^ T >] 


(A7) 


Therefore, separating the real and imaginary parts, 


f 0 (X,y)f 0 (X,ri) + f 1 (X,y)f 1 (X,Ti) = 

■^1 - cos Xa y ') (l - cos Xa^') cos >$(|y| - 111 ) (a 8 ) 


M 2 


« 2 x^ 


-f 1 (X,y)f 0 (X,r]) + f 0 ( 4 ,y)f 1 (X,Ti) = 


(ao'l 

* 2 X 2 


(l - cos Xa y ')(l - cos Xa^' ) sin X 0 (|y| - | t| | ) (A 9 ) 
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APPENDIX B 
SYMBOLS 


a velocity of sound 

a Q half root chord 

a 0 ' = a 0 /m' 

a.f. half tip chord 


a t' = a t/ m ’ 


half chord at y ^a Q - v |y|) 


a ' = a 


y “y 1 


r j m* 


aspect ratio 


(- 2 -) 

\ a o + a t/ 


half span 


b* = — ^ ° A tan A 

a o' 2P 


Cd 



Cpo 

D 


wave drag coeffigient 

two-dimensional wave drag coefficient (H) 2 / m') 
lift coefficient 
pressure coefficient 

two-dimensional pressure coefficient or assigned pressure 
coefficient 

drag 
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F(M) 
f(M) = 

G(X,q) 

g(M) = 

Ho U) ( ) 
h 0 (2) ( ) 

) 

H! (2) ( ) 
H 1 (2)( ) 

J 0 ( ) 

Jl( ) 

M 


distribution function of the harmonic source (complex); 
see equation (5) 

- J F( ) 

doublet distribution function (complex); see equation (7) 

i\G( ) 

8m' 

Hankel function of the first kind of zero order 
Hankel function of the second kind of zero order 

Hankel function of the first kind of first order 
Hankel function of the second kind of first order 

conjugate of ) (jp( ) + iY 1 ( )) 

Bessel function of zero order 
Bessel function of first order 
Mach number 


m’ = ^ M 2 - 1 

n unit vector normal to solid surface 

p pressure 

p o free-stream static pressure 


q 


velocity vector 
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r 


= \[( y - n ) 2 + z 2 


t longitudinal coordinate in transformed plane (equivalent 

to time in supersonic case) (x/m' ) 

t Q = t + a 0 * 
ti = t - a 0 ' 


u 

free-stream velocity 


u 

velocity component in 

x- (or t-) direction 

V 

velocity component in 

y-direction 

w 

velocity component in 

z-direction 

X 

longitudinal axis 


y 

spanwise axis 



y' = y/a 0 ' 


Y 0 ( ) 

Yl( ) 


a 


a. 


a_ 






Bessel function of the second kind of zero order 
Bessel function of the second kind of first order 
vertical axis 
angle of attack 

angle of attack of upper surface of airfoil 
angle of attack of lower surface of airfoil 
slope of upper surface on symmetrical airfoil 
angle of attack of cambered wing 


tan A 
3 = 

m 1 


tan A* 


* 


Po = 3 + V 
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Pi = 3 - v' 


A 

A' 


A 1 

£ 


n 

e = 

\ 


p = 


v = 


angle between midchord and y-axis (x,y-plane) 
angle between midchord and y-axis (t,y-plane) 
leading-edge sweepback angle (x,y-plane) 
trailing-edge sweepback angle (x,y-plane) 
source location along z-axis 
source location along y-axis 

tan ' 1 ?LzJL 
z 

oscillation frequency of continuous spectrum 

b/a Q 

a o " a t 
b 


I 



P 

Po 


1 - a 
b' 

source location along x-axis (t-axis) 
density 

free-stream density 


a taper ratio ^at/a Q ^ 

t = |/m' 


<t> 

X 


velocity potential 
potential in Trefftz plane 


o> = p tan A 
(% = p tan Aq 
(D i = p tan A^ 
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APPENDIX C 

INTEGRALS FOR USE IN TEXT 


The following integrals are those which are used often in the text 
and many of which are not available in the ordinary handbook. In the 
first 13 formulas, the upper line corresponds to 3 < 1, and the lower 
line to 3 > 1. 


f 

I du sin- 1 + p) = u sin" 1 + p) 


du sin" 1 (p $ |) * u sin 1 (p + ~) 


±K 


[(id- e) 2 £] 

(l - p2)V2 

cosh*” 

+K 

n ■( *1 “1 

} T (p2 - ^ 

(p 2 - l) 1 / 2 

bin 

+K 

cosh-1 

[(* - * »] 

(i - e 2 ) 1/2 


+K 

fiT « “1 


( g 2 . !)l/2 




du cosh" 1 + p) = u cosh- 1 + p) 


“—r — 7— sin * (J + fl - a2)u 

1 - p2)l/2 L 1 P kJ 


(p 2 


fjUS C06h ‘ 1 [ P T (p2 - i} k1 


du cosh -1 (p + = u cosh” 1 ^P + ~ j + < 


vtv * '“‘'i'* 8 -■>»"] 
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±K 2 [ 

[p i (i - P 2 )|]\Jl - [p ± (l - P 2 )^ 2 + sin-1 [p ± (l - P 2 )£|j 


2(1 - p 2 ) 3 / 2 1 

du \|(K + pu) 2 - u 2 = A 

f 



TK 2 1 

[p * (P 2 - i)|\J[p ♦ (P 2 - 1)^] 2 - 1 - cosh-1 |; f („ a . x) u] 

L 

2(p 2 - l) 3 / 2 | 


”1 

r 

K 2 

((1 - p 2 ) 1 1 3] \ji 


2(l - P 2 ) 3 ^ 2 

du ^|(pu T K) 2 - u 2 = i 


r I r 


K 2 

[(P 2 - l)H ! pj \j 

l 

2(p 2 - l) 3 / 2 


(e 2 -i)| + -e 


1 - cosh 


-‘f' 2 -*>!-- 


" *■[(“• e ) co8t ’' 1 (5 * p) - • o)‘ 


( K ? t) u cosh ' 1 (* * p) 


r sin-1 r P * (1 - P 2 )^l + — - [l - fp i (l - B 2 )!±l 

(1 - p 2 ) 1/2 L K -J 2(1 - p 2 ) 1 / 2 1 L K J 


2(P 2 


cosh ' 1 f T (p2 ■ l} i + g(p 2 " Ku l) i/ g J[ p ’ {p2 - - 


Jo(at)(l - cos bt) 


dt = 


kosh 


-1 b 


(a >b) 
(b > a) 


^ Jl(at) sin bt = 


UO 


b - \jb^ - 


(a > b) 
(b > a) 


Jl(at)(l - cos bt) = 


2 cosh 




(a > b) 
(b > a) 


■ Yo(at) sin bt - 


UO 


(a > b) 
(b > a) 
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^ Y 0 (at)(l - cos bt) = 


UO 


K-- 1 i) 


Q Y x (at) cos bt = < 


JO 


Vb27a2 


CM 

w 

(a £ b) 

[ cosh_1 (S] 2 j 

(b £ a) 

(a > b) 


(b > a) 



Yi(at) sin bt 
t 2 


UO 


i-^b \J b 2 - a 2 - a 2 cosh” 1 ^ 


(a > b) 


(b > a) 


r 


i 2s (a 2 - z 2 ) l/2 


(2s - l)’.a 2s+g „ 

4 S ( s + i):(s - l): 2 


^ Y 1 (at)(l - cos bt) = 


0 (a > b) 

- - a 2 

■ ~ (* > a) 




\J(t x - 3 x ti) 2 - (y - n) 2 
(y - t)) 2 


§ *, cosh-i 

v 2 v y - n 


+ < 


gsTl cosh -i (31 2 - 1)^ y - g ltl 
' h - f>iy 

sin-! K - ^ + y - 

h - Piy 


dJi 


^(t x - 3!n) 2 - (y + n) 2 
(y + n) 2 


fir 

V* 2 


= - 


fir , *1 " M 

, — - 3-i cosh”- 1 - 

v' 1 y + t) 


<?— . .O.t-1 K - - l y * ^0 

1 t x + 3 x y 


S i n -i (^-^-(^^ 1 ) 
ta+.Piy 


In the last two formulas, use upper line if pi > 1, lower line if ^ < 1. 
Similar integrals for p 0 are obtained by substituting t Q for ti, 

Po for Pi* 
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TABLE I 

INTEGRATION LIMITS OP 1) AND CORRESPONDING VALUES OF THE INTEGRAL A 1 
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TABLE II 


INTEGRATION LIMITS OF 1) AND CORRESPONDING VALDES OF THE INTEGRAL C 1 


Po > 1 (The integral exists only when t Q > y; the lower limit is. always zero.) 


Condition 


Upper 

limit 


Integral C 


to - y 


1 + Po 


< b 


to - y 


i + Po 


^ tQ ~ ^ + p 0 cosh -1 — - \|p 0 2 - 1 cosh -1 
y y 1 


-(y + Poto) 


t 0 + P Q y 


~ ~ + p 0 cosh -1 — - \|p 0 2 - 1 cosh -1 

v V * 


-(y + Poto) 


t 0 + P 0 y 


t 0 - y >, 


1 + Po 


\/(to - 3o*) 2 - (y + *) 2 „ t 0 - Pob iJ r^- Tr (Po 2 - l)b - (y + P 0 t 0 ) 

— - Po cosh 1 UPo^ - 1 cosh 1 - - 

y + b y + b 1 ° t 0 + P 0 y 


P 0 < 1 (The integral exists only when t > y; the lower limit is always zero.) 


Condition 


Upper 

limit 


Integral C 


t Q - y 


1 + Po 


< b 


t 0 - y 

1 + Po 


§ \JTT77 + £HZ + p C08h -1 h _ Bln -1 

2 N y 0 y 'f t 0 + p 0 y 


£ZIE + Po cosh-1 *£ - ^TTi7 Bln -1 L!J 

y y 1 t n + 


Poto 


P 0 y 


-^irh 


1 + Po 


J(to - Po*>) 2 - (y + t>) 2 1 to - P 0 b r — 72 , l 13 ! 1- to 2 ) - (y + Poto) 

I 1 i Po cosh -1 + \jl - Po d sin -1 - - 

y + b y + b ’ ■ t 0 + p 0 y 


table for integral D may be obtained from table II by substituting pj for p Q and 


ti for to* 
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U(M > 1) 



x = m't 
| = m'T 

(a) In physical plane (x,y). 

Figure 1.- Sweptback wing in physical and transformed planes. 

PARAMETER RELATIONS BETWEEN PHYSICAL PLANE (x,y) 

AND TRANSFORMED PLANE (t,y) 


a half root chord 
o 

aj. half tip chord 
b half span 


a t' = at/ 1 


1 + a 


A tan A 


tan a 


tan A 


1 - a _ 2p l-o 
b ' A tan a 1 + a 


ay = a Q - v|y| 


i- = ay/m' 


a '= B^/ a 


= 1 - v»b 1 


*0 + ^ 


m* (1 + a) 


m 1 = V M* - 1 


= p - v 1 


a Q surface slope 


« = p + v 




(b) In transformed plane (t,y). Mach line at 45 
Figure 1.- Concluded. 
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Figure 2.- Wave drag coefficients for the family of sweptback wings with taper 

ratio of 0.2. 
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Figure 5.- Location of cusps (p Q - 1 and ^ = 1) of the wave drag parameter in terms of 

and A tan A at o = 0.2 and 0.5. 




(a) Trailing edge at 15° from leading edge which is normal to direction 

of flight. 

Figure 7.- Downwash distribution for three infinite half -wings. 






y/t 


(b) Trailing edge at 30° from leading edge which is normal to direction 

of flight. 

Figure 7.- Continued. 


(c) Trailing edge at 45 u from leading edge which is normal to direction 

of flight. 

Figure 7.- Concluded. 
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Figure 8.- Downwash distribution over a wing tip of unit chord with a raked 

angle of 30° for three values of t. 
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Figure 9.- Downwash distribution near the nose of an infinite constant -chord sweptback wing. 
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Figure 10.- Downwash distribution at different values of y of a finite triangular wing with a 60° sweptback 

angle. 
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Figure 11.- Downwash distribution at different values of y of a tapered sweptback wing. 
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Figure 12.- Downwash distribution at different values of y of a tapered sweptback wing far downstream. 



